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Abstract: Measurements of the drift velocities of electrons and holes as functions of electric 
field strength and temperature in high-purity n- and p-type silicon with <100> crystal orientation 
are presented. The measurements cover electric field strengths between 2.4 and 50 kV/cm and tem¬ 
peratures between 233 and 333 K. Two methods have been used for extracting the drift velocities 
from current transient measurements: A time-of-flight (tof) method and fits of simulated transients 
to the measured transients, with the parameters describing the field and temperature dependence of 
the electron and hole mobilities as free parameters. A new mobility parametrization, which also 
provides a better description of existing data than previous ones, allowed an extension of the clas¬ 
sical tof method to the situation of non-uniform field sfrengfhs. For the fit method, the use of the 
convolution theorem of Fourier transforms enabled us to precisely determine the electronics trans¬ 
fer function of the complete set-up, including the sensor properties. The agreement between the 
tof and the fit method is about 1 %, which corresponds to a time-of-flight uncertainty of 30 ps for a 
pad diode of 200 pm thickness at the highest voltages. Combining our results with published data 
of low-field mobilifies, we derive parameterizafions of the drift velocities in high-ohmic <100> sil¬ 
icon for electrons and holes for field sfrengths between 0 and 50 kV/cm and temperatures between 
233 and 333 K. 
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1. Introduction 

The understanding and accurate simulation of silicon sensors requires the precise knowledge of 
basic material parameters, like the drift velocity of the charge carriers as a function of the electric 
field strength. At present, most silicon sensors used for tracking in high energy physics or for 
imaging in X-ray science use silicon with <100> lattice orientation, for which the field dependence 
of the drift velocity is only poorly known HI], 0], and parameterizations based on measurements 
exist only up to 20 kV/cm. Therefore, we performed a precision measurement of the drift velocities 
in <100> silicon for electric field strengfhs befween 2.4 and 50kV/cm and temperatures between 
233 and 333 K, which covers the operating conditions at high luminosity colliders. This paper 
summarizes the main results of our study. More details are given in Refs. 0,0,0. 


2. Mobility Parametrization 


As we have not been able to describe the measurements, which are reported in the next section, 
with the commonly used Caughey-Thomas (CT) parametrization [0, we investigated the time-of- 
flight (tof) data published by Canali et al. [0 for both the <111> and the <100> crystal directions. 
We found that a new parametrization of the carrier mobility jU(£') as a function of the electric field 
sfrength E with a threshold field £"0 improves the description of the data: 




l/Aijf* E <Eq 

+ £>£o, 


(2.1) 


with the low-field mobility /To and the saturation velocity Vsat- The root-mean-square deviation 
between fit and data is reduced by about a factor two using Eq. PTT] 

The minimum field sfrengfh of our measurements is 2.4 kV/cm, which is above the threshold 
field sfrength £0 ~ 1-8 kV/cm for electrons at 300 K, as determined from the data of Ref. [0]. 
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Thus the parametrization of the mobility used for electrons is identical to the Trofimenkoff [0] (Tr) 
parametrization: 




( 2 . 2 ) 


For holes the description of the data is significantly improved when a quadratic term is introduced: 


i/^h{E) = 


1/Mo 


E <Eq 


l/li‘i+b-{E-Eo)+c-{E-EoY E>Eo. 


(2.3) 


The temperature dependence of the individual parameters pari of Eqs. P3| is expected to follow 

a power law: 


pari{T) = pari{T = 300K) 


r[K] 

300 K 


(2.4) 


3. Measurements and analysis methods 

Sensors The sensors investigated are three pad diodes (p+-n-n+ and n+-p-p+ ) produced by two 
vendors on <100> silicon with bulk doping of 0.8 to 3.6T0^^ cm^^. Their thicknesses are 200 and 
287 p,m, which were determined to an accuracy of 1 %. Their areas are 4.4 and 25 mm^. A detailed 
description of the individual diodes is given in Ref. [0. 

Measurement set-ups The transient current technique (TCT) was employed to measure the 
charge carrier drift velocities. Sub-nanosecond laser light pulses were used to generate electron- 
hole pairs in the sensors, which drift in the electric field to the electrodes. The current induced 
by the drifting charges is amplified and recorded by an oscilloscope with 5 Gsamples/s. Light 
with a wavelength of 675 nm, which has an absorption length of a few pm in silicon, has been 
injected from both sides in order to produce signals dominated by either electrons or holes. Light 
with a wavelength of 1063 nm, which has an absorption length larger than the sensor thicknesses, 
produces signals with approximately equal contributions from electrons and holes. Most measure¬ 
ments were performed at the University of Hamburg, and control measurements were also made at 
the CERN-SSD Lab [0. 

Fit method In the fit method, the current transients are simulated with the mobility parame- 
terizations discussed in Sect. 0, convolved with the electronic transfer function and fitted to the 
measurements for determining the values of the mobility parameters. Electron-hole pairs are gen¬ 
erated according to the temperature-dependent light attenuation on a 50 nm grid, the time step for 
the drift is 10 ps, and diffusion is taken into account using the Einstein relation. A uniform doping 
is assumed resulting in an electric field with a constant gradient normal to the surface the pad diode. 
Eigure |T] shows the comparison of the fit results to the measured transients for a 200 pm thick n- 
type diode using the 675 nm laser for bias voltages between 100 and lOOOV and temperatures of 
233 and 333 K. The agreement is excellent with the exception of the hole drift data at 100 V, which 
is just above the depletion voltage of this diode. Here plasma effects llTO] become relevant due to 
the high carrier density in very low field. 
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(a) (b) 



(c) (d) 


Figure 1: Measured current transients (crosses) and the simulated ones (solid lines) using the 
parameters of Tab. |I](a) for a 200 |J.m thick n-type diode using 675 nm laser light and bias voltages 
between 100 and 1000 V. (a) Electron signal at 233 K (front-side illumination), (b) hole signal at 
233 K (rear-side illumination), (c) electron signal at 333 K, and (d) hole signal at 333 K. 


Time-of-flight method A time-of-llight (tof) method is used to cross check the fit method. The 
measured transients with the 675 nm laser light are interpolated and the difference between the 
time of the maximum slope at the rise and at the fall of the pulse, is determined. For the mobility 
parametrization l//i(£') oc a\+a 2 -E, as in Eqs. pT| , the electric field strength corresponding 
to the drift velocity v{Etof) = ^ is l/Etof = {l/E{x)), where w denotes the sensor thickness. 
For uniform bulk doping, where the electric field sfrengfh depends linearly on fhe disfance fo fhe 
elecfrodes, Etof = {Emax — Emin)/^'^{Emax/Emin), wifh Emax the maximum and Emin the minimum 
elecfric field sfrengfh in fhe sensor. 

4. Results 

The fof and fif resulfs agree wifhin 1 % for bias volfages approximafely 30 V above fhe depletion 
volfage. The resulfs obfained from fhe differenf sensors and fhe differenf sef-ups agree wifhin 4 %. 
The agreemenf befween fhe resulfs from fhe global fif for all femperafures simulfaneously fo fhe 
fifs af fhe individual femperafures is wifhin 2 %. Table |T](a) gives fhe paramefers from fhe global 


- 3 - 






































































































E = (2.4-50)kV/cm 

pari 

pari{T = 300 K) 

a, 

Md 

1430cm^/Vs 

-1.99 

Vlat 

1.05 • lO’cm/s 

-0.302 

Mo 

457 cm^/Vs 

-2.80 

b 

9.57 • 10^* s/cm 

-0.155 

c 

-3.24-10-13 s/V 

- 

Eo 

2970 V/cm 

5.63 


E = (0-50)kV/cm 

pari 

par i{T = 300K) 

at 

Md,7ac 

1530cm3/Vs 

-2.42 

^sat 

1.03 • lO’cm/s 

-0.226 

^^0Jac 

464cm3/Vs 

-2.20 

b 

9.57 • 10-^ s/cm 

-0.101 

c 

-3.31 -lO-i^s/V 

- 

Eo 

2640 V/cm 

0.526 


(a) 


(b) 


Table 1: Mobility parameters for < 100> silicon obtained from the fit to the data using Eq. for the 
electron mobility, Eq. for the hole mobility, and Eq. for the temperature dependence. The 
mean electric field sfrengfh in the sensors was (2.4 — 50) kV/cm and the temperature (233 — 333) K. 
(a) Parameters from the global fit of the data, (b) Parameters from the global fit with and the 
corresponding a, of Ref. OH] fixed. We recommend fo use (a) for E > 2.4 kV/cm. 


fif. We esfimafe thaf for electric field sfrengfhs befween 2.4 and 50 kV/cm fhe uncerfainfy of fhe 
drifl velocifies calculafed using fhese paramefers is 2.5 % for elecfrons and 5 % for holes. The 
comparison of fhe measured fransienfs fo fhe simulated ones using fhe paramefers from fhe global 
fif for fhe lowesf and highesf femperafures for differenl bias volfages is shown in Eigure |I]. 

In order fo provide a paramefrization for field sfrengfhs befween 0 and 50 kV/cm we performed 
a fif in which the low-field mobilities fto’yac *^beir temperature dependencies a, have been fixed 
fo fhe values given in Ref. [HH. The resulfs are given in Table |I](b). As there is some tension 
between the low-field mobilifies from Ref. m and our values, we recommend fo use fhe values 
given in Table |I](a) for E > 2.4 kV/cm. 

Eigure 0 compares the drift velocities calculated using the parameters of Table |I](a) for the 
<100> direction to the values of Jacoboni et al. HH] for the <111>, of Becker et al. [0] for the 
<100>, and the tof data of Canali et al. H?]] for the <100> direction. As the differences between our 
<100> results and the <111> results are typically 10 % for electrons and up to 20 % for holes at 
high electric field sfrengfhs, fhey should nof be used for simulafing <100> silicon. However, this is 
usually done, because for strong fields no parameterizafions for the <100> direction are available. 
The <100> values of Becker et al. agree with our results within 6 %, except for the hole velocity 
at strong fields, where fhese dafa had no sensitivify. The Canali fof values agree wifh our resulfs fo 
wifhin beffer fhan 10 % for E > 2.4 kV/cm excepf for fhe hole velocify al slrong fields. 

5. Summary 

The driff velocities of electrons and holes for high ohmic silicon with <100> orientation have been 
measured for electric field sfrengfhs befween 2.4 and 50 kV/cm and femperafures befween 233 
and 333 K. Compared fo previous measuremenls, significant differences of up to 20 % have been 
found. Parameterizations of the drift velocities with an estimated uncertainty of 2.5 % for electrons 
and 5 % for holes in the field sfrengfh and femperalure range of fhe measuremenls are provided. 
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Figure 2: Comparison of the drift velocities of electrons (e) and holes (h) from the global fit results 
of the data to the published data. The upper part shows the drift velocities and the lower part their 
ratios for (a) 245 K, and (b) 300 K. 


In addition, a less precise parametrization for electric field strengths between 0 and 50 kV/cm is 
presented. 
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